Distribution of maximum velocities in avalanches near the depinning transition 
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We report exact predictions for universal scaling exponents and scaling functions associated with 
the distribution of the maximum collective avalanche propagation velocities Vm in the mean field 
theory of the interface depinning transition. We derive the extreme value distribution P{vm\T) 
for the maximum velocities in avalanches of fixed duration T, and verify the results by numerical 
simulation near the critical point. We find that the tail of the distribution of maximum velocity for 
an arbitrary avalanche duration, Vm, scales as P{Vm) ~ for large These results account for 
the observed power-law distribution of the maximum amplitudes in acoustic emission experiments 
of crystal plasticity, and are also broadly applicable to other systems in the mean-field interface 
depinning universality class, ranging from magnets to earthquakes. 
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Avalanche phenomena have been observed in a wide 
variety of disordered systems that exhibit crackling 
noise near a depinning transition. Examples include 
Barkhausen noise in soft magnetic materials (TJ [2] , elas- 
tic depinning of charge density waves O H] , dynamics of 
superconductors [S], seismic activity in earthquakes [B], 
acoustic emission in mesoscopic crystal plasticity [7], 
and fracture propagation [8j. Although these varying 
materials/systems have different microscopic details, on 
long length scales, the statistical scaling behavior of 
avalanches appears to be universal. For example, the 
distributions of avalanche sizes in sheared crystals and 
in slowly magnetized soft magnets are both captured by 
the mean field theory of a slowly driven elastic interface 
in a disordered medium [U [6l [9Hl2] . 

Recent experimental studies of slip avalanches in meso- 
copic crystal plasticity have reported that the distribu- 
tion of the maximum amplitude Am of the acoustic emis- 
sion (AE) signal from each avalanche follows a power law 
P{Am) ~ ^m'", where the exponent /.t sa 2 ITBHTT] . 
Since each avalanche contributes with only one maximum 
amplitude to the histogram, many events are required to 
obtain good statistics for P{Am)- Thus the variations in 
the experimental values of /i depend on the experimental 
statistics. Owing to the proportionality between the AE 
amplitude Am and the collective velocity w,„ of disloca- 
tions [15, the distributions P{Am) and P{vm) should be 
characterized by the same scaling exponents and scaling 
functions. So far, a theoretical prediction for the value 
of the exponent ^ has been lacking. 

In this Letter, we present the first theoretical calcula- 
tion of the maximum velocity distribution, establishing a 
connection to the known classes of extreme value statis- 
tics of correlated variables. In particular, we derive the 
distribution of maximum velocities P(vm) from a mean 
field interface depinning model. We first show that the 
probability distribution function (PDF) of the maximum 
velocity for avalanches of fixed duration T follows a uni- 



versal scaling form P(w,„|r) = {2vmT)~^^^F{y/2vm/T), 
with a scaling function F{x) that can be derived exactly 
by a mapping to an equivalent problem of random excur- 
sions of Brownian motion in a logarithmic potential. Al- 
though a general theory of extremal statistics for strongly 
correlated variables is not known, much progress has been 
made already for several classes of power-law correlated 
noise with an l/w" (where w is the frequency) power 
spectrum. Brownian noise corresponds to the particular 
case where a = 2 [18l [19] . The extreme value statistics of 
power-law correlated noise typically have a robust scal- 
ing form, but the scaling function depends on boundary 
conditions, the value from which the maximum is mea- 
sured, as well as other constraints on the time evolution. 
For example, different scaling functions are obtained for 
the maximum heights of periodic Gaussian interfaces: if 
the maximum is measured relative to the spatially aver- 
aged height, the corresponding EVS is determined by the 
so-called Airy distribution function [TSJ , whereas 

measuring the maximum relative to the boundary value 
leads to the Rayleigh distribution [19l [23]. Here we 
demonstrate that our problem of maximum heights of 
amplitudes of mean field avalanches is equivalent to a re- 
lated problem whose exact solution obeys the same scal- 
ing form with a distinct function. Finally, we show that 
the overall distribution scales like P{vm) ^ by inte- 
grating P{vm\T) against the duration PDF FiT). The 
results of this study are expected to be broadly applicable 
to plasticity, earthquakes, Barkhausen noise in soft mag- 
nets and many other systems in the mean field interface 
depinning universality class. 

Our starting point is a zero-dimensional model of a 
slowly driven elastic interface in a disordered medium, 
also known as the Alessandro-Beatrice-Bertotti-Montorsi 
( ABBM) model [21] , corresponding to the dynamics of a 
particle pulled by an elastic spring and an external field 
through a random force landscape. The position of the 
particle u{t) corresponds to the center-of-mass displace- 
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ment of the interface u{t) = L~'^ J df^x u{x,t), given the 
local displacement u{x, t) at position x along the inter- 
face length L and time t for an interface of dimension d 
embedded in a (d+ l)-dimcnsional space. In the ABBM 
model, the evolution of the particle velocity v — du/dt is 
obtained by a time-differentiation of the mean field equa- 
tion of motion of the interface and given as [101 [Ml HI] 



dv 
dt 



-kv - 



(1) 



where c is the constant drift velocity, k is the elastic 
coupling constant and ^(t) is Gaussian white noise with 
autocorrelation {£,{t)(,{t')) = 2DS{t — t') where _D is a 
constant measure of disorder. In numerous studies, this 
model has been shown to reproduce well the universal 
scaling laws for the size and duration distributions near 
quasi-static depinning for systems with long-range inter- 
actions [U Uni [21H1B] . It is relevant for the calculations 
below to recall the power law decay of the distribution 
of avalanche durations Ft{T) « r-(2-<=/-D)/j, (^), with 
a rate-dependent exponent [TJ[M1[5S], and such that, for 
T <^T* and c — > 0, the distribution of durations follows 
the mean field scaling law Ft{T) « . The exponen- 
tial scahng function Jt[T/T*) and the cutoff T*{k) to 
the scale invariance can be computed analytically in the 
limit of c = [15J [27] . The probability distribution for 
the velocity follows the Fokker-Planck equation 



dtP{v,t) = {{kv - c)Piv,t))+Ddl {vP{v,t)) 
which has a steady state solution given by 
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where T{z) is the Gamma function, c = c/D and k = 
k/D [H [Ml US [17] . The power-law exponent depends 
linearly on the driving rate c, such that in the adiabatic 
limit c — >■ 0, the distribution approaches the well-known 
scaling, which has been verified by experiments on 
the dynamics of domain walls in ferromagnets |28j . 

Maximum velocity distribution for avalanches of fixed 
duration:- By a change of variables to x = 2y/v, 
Eq. ([ij transforms to an additive-noise Langevin equa- 
tion dx/dt ^ -kx/2 + (2c - D)/x + ^{t). The addi- 
tional 1 /x term comes from the Ito interpretation of the 
multiplicative noise in Eq. ([T]). This choice yields the 
correct Eq. Thus, in the adiabatic limit, near de- 
pinning, where c — >■ and fc — >■ 0, the velocity evolution 
can be mapped onto a one-dimensional (ID) Brownian 
motion in a logarithmic potential. An avalanche of dura- 
tion T corresponds to an excursion, i.e. a path x{t) with 
x(0) = x(T) = and x(t) > for < t < T. The ex- 
treme displacement distribution for Brownian excursions 
can be derived using the path integral formalism found 
in Refs. |1TJ[53]. We adapt this method to our problem. 
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FIG. 1. (Color online) (a) Data collapse of the PDF P(Dm|T) 
from numerical integration of Eq. |T| in the Ito interpretation 
with parameter values fc = c = and D = 1/2. Large du- 
rations (T ~ 1000) are required to obtain the scaling regime 
where Eq. Q holds. The collapse fits very well with the an- 
alytically determined F{x), which is represented by the solid 
line. The inset figure shows the P{vm\T)'s for different dura- 
tions before the rescaling. In panel (b), we show {vm\T} as a 
function of T, with the solid line representing the analytical 
solution from Eq. ||8| with L> = 1/2. 



and determine the cumulative distribution CRw{xin\T) 
of the maximum displacement during excursions for a 
Brownian motion in a logarithmic trap. The cumulative 
distribution can be defined as Ciiw{xm\T) — 



lim 

£-fO 
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(4) 



where the Lagrangian is given hy Ce — 4i5(i+x)^- 
The theta function products in the numerator indicate 
that only paths that stay positive-valued between t = 
and t = T and have a maximum distance from the ori- 
gin not greater than counted. The denomina- 
tor is a normalization factor, counting any excursion of 
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duration T without regard to its maximum value. The 
Fokker-Planck equation (Eq. ([2])) with c = A: = in 
terms of the variable x is Bessel's equation of order 1, 
thus the path integrals from Eq. Q can be written as 
the matrix elements (e| exp(— /fT)|e) of the Hamiltonian 
H = —d^ — dx/x+ l/x"^ with appropriate boundary con- 
ditions, and then expanded in terms of Bessel functions 
(details are presented in [29]). From CRw{xm\T), the 
PDF P{x„,\T) = dx^C{xJ{r) is determined. We find 
that the P{x„i\T) has the scaling form 



with scaling function 

-J oo 



n— 1 ^ ^ 
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where A„ is the n'th zero of the Bessel fimction Ji{x). 
From Eq. ([5]), it also follows that the average maxi- 
mum displacement scales with duration as T^/^, like the 
average maximum relative heights of fluctuating inter- 
faces im US]- Returning to the physical variable 
w, we find that the maximum velocity distribution in 
avalanches of fixed duration T has the scaling form 




(7) 



The average maximum velocity dependence on avalanche 
duration T can be obtained as the first moment of the 
conditional distribution 



/•oo 

{v,n\T) = / dVmVmPiv,n\T) = DT , 
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where we used the fact that J dxx'^F{x) = 2. 

Using the same method, we also determine that the 
PDF of the instantaneous velocity v at time t in an 
avalanche of duration T is given by 



P{v,t\T) 



T 



-vT/(Dt{T-t)) 



(9) 



The first moment of P{v,t\T) gives {v{t)\T) = 2Dt{T ~ 
t)/T^ which is the parabolic average avalanche shape dis- 
cussed in "2", "30", "3T] . 

We have verified Eqs. ([7| and ([s]) numerically by inte- 
grating Eq (see Fig ([if). Large durations (T - 1000) 
must be explored for the scaling function to converge to 
the one predicted by our continuum derivation, but the 
results are in accord with predictions. We obtained im- 
proved statistical results for data collapse to the same 
scaling function using the computationally-efficient dis- 
crete velocity shell model [HI [31] , which obtains its scal- 
ing regime at smaller durations. These results will be 
reported elsewhere ^29,. 
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FIG. 2. (Color online) In panel (a), the PDF P(vm) obtained 
from numerical integration of Eq. ([l]) is shown for several val- 
ues of c, all with k = 10~* and D = 1/2. The PDFs are 
offset vertically so they can be clearly distinguished. Above 
each curve, a guideline is drawn indicating the power law an- 
alytically predicted from Eq. (111. In panel (b), we show the 
instantaneous velocity PDF P{v) for several values of c. The 
solid lines represent the functional form predicted by Eq. ([3| . 
Again, the PDFs are offset for visibility. 



Although our analytical calculation was performed ex- 
actly only for fc c = 0, the scaling form in Eq. ([T]) gives 
a good collapse of simulation data away from criticality 
as well, even when including durations T ^ T*{k) 20]. 
Therefore the dependence of the scaling form on k is likely 
to be weak. Since the driving rate parameter c is dimen- 
sionless, one might anticipate that nonzero values of the 
driving rate c modify the scaling function, but not the 
scaling form. Indeed, the exact form of the modification 
can be calculated analytically with a slight generalization 
of the above calculation [29] . 

Maximum and instantaneous velocity statistics:- We now 
investigate the maximum velocity distribution integrated 
over all durations. This distribution is equivalent to the 
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P{A) of the maximum AE amplitude, A, deduced from 
the time series in AE experiments on crystal plasticity. 
From Eq. ([o]), we can determine the PDF of the maxi- 
mum avalanche velocity P{vm) by integrating P{vm\T) 
over avalanche durations T weighted by their distribution 
Ft{T) - r-2+£^ foj. T <^T*. Our numerics indicate that 
P{vm\T) satisfies Eq. Q at least for durations T < T* . 
Thus, the distribution of Pivm) is 



PM 



dT 



P{Vm\T), 



and near depinning, where T* — >■ c», we have 



(10) 



(11) 



Similarly, we obtain the PDF P{v) for the instantaneous 
velocity v at an arbitrary time by integrating P{v,t\T) 
over the time spent in an avalanche of duration T and 
then over the distribution of durations, giving 



P{v)= ( dT j dtP{v,t\T)F{T) 
Jo Jo 



T* 



dT 



T3 



^Giv/DT) 



(12) 
(13) 



where G(a;) = du{u{l ~ u)) ^ exp{~x{u{l — u)) 



G{x) 



for a; ^ 1 and decays exponentially for 



The distributions agree with the predictions of Eqs. ( 11 ) 
and 

In addition to the exponents, it would be interesting 
to measure the predicted scaling form of the P(w,„|T) 
over fixed durations from AE experiments, in the cor- 
responding regime where the distribution of maximum 
amplitudes P{Am) ^ ^,7^ was observed [7l fT3HT7] . The 
distribution P{vm\T) was calculated exactly only at the 
depinning transition with k — c — but numerical evi- 
dence strongly suggests that an indistinguishable scaling 
form occurs away from the transition, whose dependence 
on the elastic coupling constant and the driving rate still 
needs to be studied in more detail. Finally, what happens 
beyond mean field theory remains an open question, de- 
spite the apparently good agreement of our calculations 
with available experimental data. For instance, it is un- 
clear to us whether the exponent ^ in P{vm) = is 
expressible in terms of the other avalanche exponents (t, 
avz, etc.) or if it is independent. 
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